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PREFACE 



It has become apparent during the course of instruction 
at the Naval Postgraduate School that the phenomena of heat 
transfer is of ever increasing importance in many fields of 
engineering. The attainment of higher thermal efficiencies 
ultimately becomes a problem in heat transfer. The concent 
of the heat transfer coefficient has become an accepted me- 
thod for predicting the relative ease with which heat flow 
will occur. 

In conversations with Professor E. E. Drucker, of the 
U. S. haval Postgraduate School, the authors were made aware 
of the lack of adeouate data on transfer coefficients as they 
relate to wet steam. Such lack of data was considered suf- 
ficiently important to be investigated. 

The experimental work of this thesis was conducted at 
the U. S. Naval Postgraduate School from March to June of 1953* 

The authors desire to express their appreciation to Pro- 
fessor E. E. Drucker for his extensive assistance in the pre- 
paration of this thesis, to the Research Division of the Bab- 
cock and Wilcox Company for their invaluable aid in the in- 
strumentation of the test section, and also to Chief Jones, 

USN, and Chief Wallace, USN, for their aid in the construction 
of the experimental set-up. 
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SUMMARY 



The object of the development and preparation of this 
thesis was to determine the local or point coefficients of 
heat transfer to wet steam at moderate pressures, and to note 
their manner of variation with respect to moisture content 
of the steam. To accomplish this, steam at a known condition 
was passed through a vertical steel tube of 0.50 inches in- 
side diameter which was heated electrically. 

Coefficients were determined under the following con- 
ditions: moisture content in the range from 0 to 6 percent, 
absolute pressures up to 215 psia, and a variation in mass 
velocity from 73,000 to 172,000 pounds per hour per souare 
foot . 

The effects of moisture content were found to be signifi 
cant, causing a rapid increase in the coefficient with in- 
creasing percent moisture. The local coefficients as they 
are herein presented vary from 110 BTU’ s/hr-Ft^-°F for satur- 

p 

ated steam to values approaching 65 OO BTU’s/hr-Ft -°F for 
steam with a moisture content of about 6 %. 
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INTRODUCTION 



Steam power has been long established as the primary 
means of ship propulsion in the U. S. Navy. In the advent 
of atomic power for naval vessels, it is entirely probable 
that steam will continue to be the means by which this power 
is transmitted to the main engines. 

The Navy’s interest in the improvement of the thermal 
efficiencies of boilers has led to the utilization of higher 
temperatures and pressures. It is felt that additional im- 
provement may be brought about by critical re-evaluation of 
the current concept of boiler design as it relates to" heat 
transfer. Consequently, it is important to know as much as 
possible about the mechanism of heat flow to the working 
medium as it passes through successive stages from satur- 
ated water to saturated steam. One aspect of this problem 
may be solved by experimentally determining the coefficients 
of heat transfer to steam with varying degrees of moisture 
content . 

In this thesis work, experimentation was conducted in 
an effort to determine the coefficients for steam with a 
quality of 96 to 100 percent. Although this is a fairly nar- 
row range of moisture content, the results did tend to indicate 
the trend with which the transfer coefficients varied with 
quality. Time was a limiting factor in preventing further in- 
vestigation of steam with higher moisture content. 
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Heat transfer coefficients of this kind, that is, coef- 
ficients relating to forced convection, are usually corre- 
lated with respect to the local Nusselt, Reynolds, and Prandtl 
numbers. Such correlation, however, is invariably made for 
fluids which are of a single phase only. Obviously this was 
not the case in the present investigation. In addition, a 
thorough survey of all existing literature failed to yield 
any data on values of density, viscosity, or thermal con- 
ductivity of wet steam, thus preventing the computation of 
the above dimensionless quantities. Further research to de- 
termine some of these parameters would be extremely useful. 

This work may be considered as an extension of the in- 
vestigation of heat transfer in two phase flow. It repre- 
sents some of the initial data available for steam at relatively 
low flow rates and low heat flux densities. 
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CHAPTER I 



REVIEW OF LITERATURE ON TWO PHASE FLOW 

A survey of the broad field of heat transfer indicates 
that a large quantity of data has been collected on the heat 
flow to fluids inside tubes. Unfortunately very little data 
has been published with respect to gas-liquid mixtures. 

The lack of information on temperature and velocity 
distributions has made it necessary to correlate experimental 
data on heat transfer coefficients by means of empirical 
equations in terms of certain dimensionless groups. Perhaps 
the most highly utilized equation is the Nusselt relation for 

single phase flow in forced convection through a tube (4): 

Nu=0.024 Re°- S Pr 0 *^ (1) 

It has been found by McAdams and others (4), (5) from 
data for vaporization inside tubes that the local film co- 
efficients, in the preheat section where the fluid is ini- 
tially heated, are usually somewhat higher than would be pre- 
dicted from the Nusselt relation. This result is logical 
from the standpoint that some vapor bubbles may have formed 
in the preheating section, with the possibility that they 
disrupt the film locally as they are formed on the tube wall. 
As they move out into the fluid they would condense thus re- 
leasing their heat of condensation. The combined effect is 
an increase in heat transfer coefficient. In a comparable 
manner, it might be expected that as far as the influence on 
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the film coefficient is concerned, the coefficient will in- 
crease when gas passes through a system simultaneously with 
a liouid. 

In 1949 Bergelin(2) observed that if a gas and a liquid 
pass simultaneously in upward flow through a tube, as the 
gas-liauid ratio is increased, the following three flow con- 
ditions could be distinguished: 

(a) bubble flow, when gas bubbles pass individually 
through the tube; 

(b) slug flow, when alternatively slugs of gas and 
liouid pass through the tube, and 

(c) annular flow, when the liouid flows in an annulus 
along the tube wall, while the gas passes at a much higher 
velocity through th« center of the tube. Verschoor and 
Stemerding( 8 ) in their investigation of the transfer coeffi- 
cients to air and water mixtures verified the observations 

of Bergelin. They found that as their air to water mixture 
was increased, the heat transfer coefficient increased until 
a maximum was reached that coincided with the transition from 
slug to annular flow. Prior to reaching this maximum there 
was a change in slope in their curve which corresponded to 
the change from bubble flow to slug flow, it is interest- 
ing to note that their heat transfer coefficients increased 
with an increase in flow density, however this was of second- 
ary importance compared to the influence of varying the air- 
water mixture. 
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Yoder and Dodgp(9j measured film coefficients of Freon-12 
boiling in a vertical tube. Their results showed that pres- 
sure and flow rate had very little effect on the coefficients, 
but as the relative amount of vapor increased from 4 0 to 100 
percent, there was a very significant decrease in the heat 
transfer coefficients. 
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CHAPTER II 



EQUIPMENT 

The equipment was a flow system in which steam was taken 
from a main line and passed through a throttle valve, a cooler 
the test section, a condenser, and to a weigh tank. The as- 
sembly is shown schematically in Figure 1. The steam was sup- 
plied from a Babcock and Vilcox Fii boiler at a pressure of 
200 psig. After passing through a throttle valve the steam 
was made to flow through a cooler which was intended for use 
as a method of varying the ouality. It was never subse- 
quently used since time limited the investigation to a per- 
cent moisture not greater than 6 percent. On either end of 
the test section an Ellison combination throttling and separat 
ing calorimeter was installed for quality determination. Prior 
to the condenser an orifice of 0.166 inches diameter was in- 
stalled to maintain the pressure in the test section. A by- 
pass was provided around the orifice to permit adjustment of 
the flow rate independently of pressure. 

The test section as seen in Figures 2 and 4 consisted of 
a vertical steel tube having an inside diameter of 0.50 inches 
and an outside diameter of 1.25 inches. Its overall length 
was approximately 27 g inches, but with an active length of 
24 inches for purposes of heat transfer. Recesses were cut 
in either end of the test section to minimize axial heat loss. 
The thermocouples in the test section were made of No. 28 
gauge chromel-alumel wire, and were installed by the Research 
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Division of the Babcock and Wilcox Company. Four pairs of 
thermocouples were installed along the tube, one of each pair 
at the outer surface of the tube, the other at a known depth 
below the surface thermocouples. The leads of the thermo- 
couples were threaded through a milled slot and out a i n 
tube welded to the test section. The original thought was 
that the thermocouples as located would permit determination 
of the temperature gradient across the wall of the tube, thus 
the inside surface temperature could be found by extrapola- 
tion. 

At either end of the test section, thermocouple probes 
were projected into the fluid stream for measuring the tempera- 
ture of the steam. TKese probes were obtained from Leeds and 
Northrop, and manufactured of iron-constantan wire. 

The heating element consisted of four independent heaters, 
each containing 29. # feet of No. 1$ Nichrome V wire. The 
wire was wound around the tube in a single layer. Each heater 
was approximately 6 inches in length. In order to prevent 
shorting of the heater coils, the wire was insulated with glass 
sleeving, and a thin layer of mica, 0.0015 inches thick, was 
wrapped around the tube before the heaters were wound. Power 
supply to the heaters was controlled by four 115v-15amp Power- 
stats. See Figure 3. The total power to each heater was 
measured by portable wattmeters. 

The insulation of the test section was accomplished by 
winding the tube, outside the heaters, with strips of 1/8” 



thick asbestos cloth. This was followed by applying a length 
of standard l| n magnesia pipe insulation. The entire section 
was then covered with an aluminum sheet casing. Provisions 
were made so that air, heated to approximately the same tem- 
perature as the outer surface of the insulation, could be 
supplied to the space between the insulation and casing, thus 
minimizing any heat loss due to convection, and reducing the 
conducted heat loss through the insulation. 
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CHAPTER III 



OPERATING PROCEDURE 

There were four primary variables associated with the 
apparatus: heat flux, flow rate, duality, and pressure. It 

was found, however, during preliminary experimentation, that 
the minimum moisture content of the steam available for use 
was approximately 6%. Interest at this time was primarily 
centered on the region between 0 and 10 percent moisture, 
thus necessitating throttling the steam to lower pressures 
in order to obtain a given quality at the entrance of the test 
section. Such an arrangement meant that only the heat flux, 
pressure, and flow rate could be varied independently of one 
another. The only other limiting factor associated with the 
equipment was, that due to the design of the heaters and the 
power supply used, the maximum heat input was restricted to 
1 kw per heater, or a total of 4 kw. 

The procedure was to select an approximate quality at 
the entrance of the test section by setting the throttle valve 
accordingly. The flow rate was either taken as that provided 
by the pressure across the orifice in the line to the con- 
denser or adjusted by means of the by pass around the orifice. 
Power was then supplied to the heaters and the system was al- 
lowed to come to a steady-state condition. 

The approach to steady-state condition was determined 
by noting the change in the millivolt reading of the thermo- 
couples contained in the test section. \Ihen they became 
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approximately constant, steady-state was considered to exist. 
It normally took from 1/2 to 3 /K hours before the system 
reached equilibrium. It was never possible to keep tempera- 
tures constant over a period longer than a few minutes, even 
though heat input and flow rate appeared unchanged. Each 
run was made at constant flow rate, heat flux, pressure and 
quality of the inlet steam. 

For each run, temperatures, pressures, flow rate, and 
power input were recorded. Data points were taken at various 
settings of power inputs, flow rates, and pressures. Most 
of the runs were concluded at a power input such that the 
minimum change in quality which occurred through the test 
section was yjo. 
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CHAPTER IV 



METHOD OF CALCULATION 

The fundamental law of heat transfer by convection may 
be stated by the mathematical expression: ( 4 ) 

p^/?/74r ( 2 ) 

The factor of proportionality h, defined in the equation, is 
called "the film coefficient of heat transfer." It is this 
coefficient which was computed. 

For each set of conditions imposed on the test section, 
the density of heat flux, the inside surface temperature of 
the tube, the ^3 T, the quality of the steam, and the local 
coefficient of heat transfer were calculated. Each of these 
factors was evaluated at four positions corresponding to the 
locations of the thermocouples inside the tube wall. 

The power dissipated in the heating elements was ob- 
tained from wattmeter measurements. The heat flux density 
was then calculated from the power dissipated and the dimen- 
sions of the heat transfer surface. 

Early in the process of collecting data for the investi- 
gation, it was discovered that the thermocouple pairs in the 
test section were indicating unreasonable temperature gradi- 
ents through the tube wall. It became necessary to calcu- 
late the inner surface temperature by applying to the depth 
thermocouple a computed temperature drop corresponding to 
the heat input. The temperature drop was calculated using 
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the Fourier conduction equation for steady-state heat trans- 
mission: (4) 

£ = ~ A/> $r ( 3 ) 

For a circular tube, this eciuation may be integrated to give: 

O = 2TT/rL AT (4) 

The depth thermocouples were shielded from any direct influ- 
ence of the heaters, and thus were considered to be more ac- 
curate than those on the surface of the tube. 

The temperature of the steam was measured at both the 
entrance and exit of the test section, however, these measure- 
ments were consistently erratic and failed to give desired 
results. It was felt that greater reliability would result 
from assuming the saturation temperature corresponding to the 
absolute pressure in the test section for computing the ^1T. 

In the preliminary phase of the investigation it was 
discovered that the reproducibility of the calorimeters, for 
determination of quality below the level which could be found 
by the throttling principle alone, was rather poor. With re- 
sults taken over a long period of time it was found that the 
average quality of the steam at the test section, at a pres- 
sure of 200 psig, was 94$ • Further use of the calorimeters 

was discontinued, and this value of 94$ for quality was used 

\ 

as a basis for all further computations. The quality of the 
steam, at each point in the test section corresponding to the 
location of the thermocouples, was computed by adjusting the 
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enthalpy of the incoming steam in an amount equal to the 
heat input to the test section. The quality at pressures 
lower than 200 psig was determined by the use of the Ilollier 
Diagram. 
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CHAPTER V 



CONCLUSIONS AND RESULTS 

One hundred runs were made of which forty-six were suc- 
cessful. The remainder of the runs were eliminated because 
the tube wall temperatures obtained indicated that steady- 
state conditions had not yet been reached. The ranges of 
pressure, temperature, quality, and mass rate of flow are 
given in Table I of Appendix I . 

In the process of collecting data a number of experi- 
mental errors became apparent. These errors, however, did 
not exceed the range of expectancy ordinarily associated with 
this type of investigation. The temperature determinations 
were expected to be in error at least 2 to 3 %. This was not 
serious as long as the temperature difference between the 
fluid and the wall surface was large. At small temperature 
differences care had to be taken to insure that unreasonable 
results were not obtained. Another error is introduced in 
the calculation of the transfer coefficient by neglecting 
axial heat loss from the test section. In estimating this 
loss it was found to be small, in the range of 10 to 40 BTU/ 
hr. The error due to the radial heat loss was approximately 
2 ‘ % of the input power. 

Curves for the heat transfer coefficient versus the per- 
cent moisture were made for each mass velocity. In addition 
a plot was made of the transfer coefficient versus the tem- 
perature difference between the wall and the fluid, and one 
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of the transfer coefficient versus the mass rate of flow. 
These curves are seen in Figures 6 through 12. The physical 
properties of steam-water mixtures were not available, thus 
preventing any correlation of the transfer coefficient with 
any existing data. 

The curves depicting the heat transfer coefficient 
versus the percent moisture show a very gradual rise in the 
coefficient up to moisture content between one and two per- 
cent. Above this value the coefficient increased rapidly as 
the moisture increased from two to five percent. A possible 
explanation of this phenomenon is that at high qualities 
there are only small droplets of water in the fluid, but not 
enough to materially disrupt the boundary layer between the 
steam and the wall of the test section. As the amount of 
moisture in the steam increases, the water droplets become 
larger and begin to rupture the boundary layer and collect 
on the wall of the test section, with a consequent increase 
in the heat transfer coefficient. The sharp change in slope 
of the curve is the point at which the water begins to col- 
lect in droplets on the tube wall. With a further increase 
in moisture, there is an approach to annular flow where the 
boundary layer becomes completely liauid. At this point 
the transfer coefficient should reach a maximum and then 
start to decrease. Fully developed annular flow was never 
reached in this experiment. 
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Figure 11 shows the influence of mass flow rate on the 
heat transfer coefficient. It is surprising to note that the 
coefficient decreased as the flow rate increased to a value 
of about 200 lbs/hr. With flow rates above this value, the 
coefficient increased. Only one point on the curve was es- 
tablished beyond this transition region. However, the pre- 
dicted curve is extrapolated (dotted line) on the plot. The 
initial decrease in the coefficient was possibly due to the 
water droplets on the wall being flattened out as the flow 
rate increased and thereby forming a liquid boundary layer. 

As the mass flow rate increased beyond 200 lbs/hr the veloc- 
ity of the gas was apparently sufficient to start wiping the 
water droplets off the wall instead of just flattening them 
out. This is ho w the continued upward trend in the curve 
was predicted. 

Figure 12 indicates curves of the transfer coefficient 
versus the difference in temperature between the tube wall 
and the fluid. This shows the correlation between the heat 
transfer coefficient and temperature difference, for a given 
range in quality of about five percent. These curves appear 
to be in the region of unstable film boiling when compared to 
the boiling curve for free conveeticn(7), but the coefficients 
determined in the investigation are necessarily higher than 
those of the boiling curve due to the influence of forced 
convection . 
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The available literature and results of this investiga- 
tion permits the conclusion that a maximum heat transfer co- 
efficient exists with respect to the moisture content of the 
mixture. This point has not been within the range of the 
subject investigation. The location of this maximum repre- 
sents an optimum for future design of certain types of heat 
transfer equipment. 

Additional research would be amply justified in order 
that the peak in this curve be definitely established. 
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TABLE III 



CALCULATED TEMPERATURE DROP FOR VaRIoUS FOYER INPUTS 



Power Input per 
Coil in Kilowatts 


Temperature Difference 
Between Surface Thermo- 
couple and Inside Surface 


Temperature Dif- 
ference Between 
Depth Couple and 
Inside Surface 


300 


11.94 


2. BO 


350 


13.92 


3.26 


400 


15.91 


3.73 


450 


17.90 


4.20 


O 

O 

or\ 


19.92 


4.66 


550 


21.9 


5.13 


600 


23.9 


5.60 


650 


25.9 


6.06 


700 


27.9 


6.53 


750 


29. B 


u.OO 


BOO 


31 . B 


u .46 


B 50 


33. 8 


7.92 


900 


35. B 


B.40 


950 


37. B 


B.B6 


1000 


39. B 


9.32 
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F45 Heat transfer coefficients of 

steam water mixtures. 
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